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Abstract 

Experiments and numerical data on the correlation length £(T) for 
large S disagree strongly with the theoretical prediction based on the 
effective field theory prescription of the magnon physics. The reason is 
that for large S, at any accessible f(T), the cut-off effects from the non- 
magnon scales become large and can not be treated by an effective field 
theory. We study these effects in a spin-wave expansion. The corrected 
prediction on £(T) connects the renormalized classical and the classical 
scaling regions smoothly and comes close to the data. 



1 On leave from the Institute of Theoretical Physics, University of Bern, Sidlerstrasse 5, 
CH-3012 Bern, Switzerland. 



There exist a number of quasi-2D antiferromagnetic compounds, including 
spin-1/2 antifcrromagnets (AFMs) La 2 Cu0 4 and Sr 2 Cu0 2 Cl2 @, spin-1 AFMs 
La 2 Ni0 4 and K 2 NiF 4 [§ and spin-5/2 AFM Rb 2 MnF 4 [§ f§ § whose magnetic 
behavior is well described by the 2D quantum Heisenberg model^] 

H = JJ2 S«+?S„ , S 2 n = S(S + l). (1) 

n,i 

The low temperature properties of this model are dominated by magnon exci- 
tations and can be described by an 0(3) invariant effective field theory 
whose leading part contains two parameters only: p s (spin stiffness) and c (spin- 
wave velocity) 

Reading = £ ^ ^ J fa x) Q^ x) ^ B? = 1 , (2) 

where = 0, 1, 2 refer to r, x\, x 2 , respectively. In their workjlO|], Chakravarty, 
Halperin and Nelson have used renormalization group to connect this effective 
theory with the d = 2 classical 0(3) non- linear er-model and, among other 
results, predicted the asymptotic, small temperature behavior of the correlation 
length £ ~ c/2irp s exp(27rp s /T). The exact mass gap of the <7-model[[ll) and 
the two-loop corrections lead finally to the asymptotic expression [[l2"| 

It has been demonstrated in a careful numerical study lljj] that the explicitly 
given terms in eq. (||) are consistent with the numerical data at very large 
correlation lengths (low temperatures) for S = 1/2. At moderate correlation 
lengths, however experiments numerical data]T3]|-p6{, series expansions 
[ p"7[ and a semiclassical model^8[ indicate a significant discrepancy which is 
increasing rapidly with S. It was subsequently realized[[l9| [l?]] that the 
quantum Heisenberg model contains a sequence of crossovers (depending on T 
and the parameters J and S), and eq. (||) is valid only in a corner called the 
region of 'renormalized classical scaling'. 

This situation is unsatisfactory since the basic feature of the quantum Heisen- 
berg model in eq. (0), namely that it can be mapped onto a simpler model, is 
valid beyond the region of renormalized classical scaling. The purpose of this 
paper is to determine the corrections to eq. (^) making it applicable at moder- 
ate correlation lengths also and for any value of S including the large S limit 
('classical scaling region'). 

We calculated the corrections to eq. @ (which are due to cut-off effects in the 
quantum Heisenberg model) in leading non-trivial order of spin-wave expansion 

2 We use ft = 1, fcg = 1 convention. 
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Table 1: The correcting factor exp(— C) as the function 7 = 2JS/T. 
and obtained 

^ = 8 2^ s eXp (-Z^ i 1 - 4^1 ) C ^-°^ < ? = IT ' (4) 

where the correction factor exp(— C(j)) is given in Table 1 for different 7 values. 
Here 7 = 2JS/T ~ A cut /T(l + 0(1/ S)), where A cut = c/o and a is the lattice 
unit. In the limit T — > 0, (5 fixed), we have 7 — > 00, and C ~ 7~ 2 , leading 
to the old result in eq. For large 5 and p s ^> T ^> A cut (classical scaling 
region), we have 7 — > and C(j) = ir/2 + In 8 + In 7 + 0(7) which gives 

1 exp(-f ) e T ,27rp. / n( A cut \ , . 

-Classical = f=~ ^ exp(— — ) 1 - 0(— — ) , (5) 

a V32 8 2Txp s T \ T J 

where p s ~ JS" 2 ~ Pdass- The fact that the prefactor exp(— ir/2)/\/32 comes 
out automatically from the corrections is a non-trivial test on our r esultfl For 
intermediate 7 values, exp(— C(j)) gives a significant factor as shown in fig. [j], 
where, for correlation lengths larger than 10, the corrected theoretical prediction 
is compared with MC and series expansion data and with a semiclassical model.^ 

We shall first discuss the steps leading to eq. (|4|), then we close the paper 
with a few remarks. 

3 Eq. ([^I) has the expected form [|2c|: for 5 — » 00, the quantum Heisenberg model goes 
over to the standard Jattice regularized 2D 0(3) c-model with bare coupling gs = JS 2 /T. 
The prefactor in e q. (|5[ ) is the ratio between the renormalization group invariant scales on the 
lattice and in the MS renormalization scheme calculated long time ago in rcf . pl| . 

4 Unfortunately, some of the most interesting large correlation length MC results at S = 1 
and 5/2 |16| are not yet published and missing from this figure. 
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2D quantum Heisenberg model vs. the classical 2D 0(3) non-linear a-model: In 
order to extend the result in eq. (^) we recall that the quantum spin model can 
be mapped onto the 2D a-model 

Aa = Y g j d2x d l e{x)d l e{x) , e 2 = 1 , i = 1, 2 (6) 

under the condition that the correlation length £ of the quantum model satisfies 
£ ^> c/T and £ ^> a, where c/T is the length scale defined by the temperature. 
In particular, one is not forced to consider the parameter region only where the 
leading effective field theory in eq. (^|) correctly describes the magnon physics. 
This is important, since eq. (|J) does not give account of the cut-off effects in 
the quantum Heisenberg model. The cut-off effects are related to the fact that 
the non-magnon length scales like cj p s are not much larger than the lattice unit 
a. Actually, for large S, c/p s ~ 1/S becomes much smaller than a. The cut-off 
effects enter the effective prescription eq. (|) on the 4-derivative level first Q and 
contribute to the 0(T 2 ) correction in eq. (|3) when T — ► with J and S fixed. In 
this sense eq. (|J) and eq. (H) become correct at sufficiently small temperatures 
for any given S. For large S, however this happens only at astronomically large 
correlation lengths, while the mapping to the 2D cr-model is valid much earlier. 

Let us turn now to the argument concerning the mapping of the 2D quantum 
Heisenberg model onto the classical 2D tr-model.^ The partition function of the 
quantum Heisenberg model can be represented in terms of a path integral using 
coherent states, for example [p3|. The action is constructed in terms of a 3- 
component classical field e(r, m, 712), e 2 = 1, where n±, n2 are coordinates of the 
two-dimensional spatial lattice and r S (0, c/T) is the continuous coordinate of 
the periodic imaginary time direction. If the correlation length £ is much larger 
than c/T (which is the case for small temperatures), in units of £ we have a 
thin slab with two infinite space directions. This is just a special regularization 
of the 2D non- linear cr-model in eq. (||) . 

Connecting the parameters of the quantum and the classical models: The map- 
ping discussed above leads to quantitative predictions if we connect the param- 
eters of the quantum Heisenberg model with the coupling g of the non-linear 
cr-model. A possibility is to choose a convenient long-distance quantity and 
calculate it in both models. The predictions should match leading to the con- 
nection we are looking for. 

A convenient low-energy quantity is the free energy density f as the function 
of the chemical potential hF\ This was the choice in ref. 1 12 , where f(h) 



/(0) was calculated in two-loop perturbation theory in the 2D classical cr-model 

5 For related arguments and different wording, see Jl2| 

6 Unlike the magnetic field, the chemical potential is rcnormalization group invariant. 
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eq. (^|) and also in the effective model eq. (||) (after introducing a chemical 
potential in both models in an equivalent way). Matching the results gives Q 

= ^+« + /3-, a = 0,/3 = ~— Lj, (7) 
9{P)\ p =t/c T Ps 4(2tt) 2 

where g(p) is the renormalized coupling constant of the er-model at momentum 
p. Using the relation p"l|] 

e 1 q(p) , 2tt , / 1 , „ \ 
e=--^exp T - 1 + - 5 W + ... , 8 

eq. (||) leads to the prediction in eq. (||). 

Calculating f(h) — /(0) in the quantum Heisenberg model with the help of 
the effective model in eq. (rather than using the Hamilton operator directly) 
is an elegant, powerful method. The underlying assumption is that the higher 
derivative terms in the effective action A c r, which do not contribute up to the 
two-loop order|9], give controllable, small corrections. As we discussed before, 
due to the large cut-off effects at large 5 in the quantum model, this assumption 
is true at very low temperatures (astronomically large correlation lengths) only. 

In order to calculate the corrections to eq. (^) due to cut-off effects, we have 
to abandon the effective prescription^] and work with the Hamilton operator 
directly. We shall use spin-wave expansion (SWE) to derive f(h) — /(0). It is 
natural to consider T an 0(5) quantity when doing thermodynamics in SWE. 
Then the 3 terms in eq. (0) are 0(5) , 0(1) and 0(1/5). Our SWE runs up to 
second order and so identifies corrections to the first two terms in eq. (0). 

The chemical potential introduced in the classical a-model (see, eq. (4) in 
rcf . Jl2[ ) corresponds to an imaginary twist in the quantum Heisenberg model^J 

H(h) = H(0) + J^{(cosh(/ l a)-l)(S I 1 l+i S J 1 l + S3 +i SfJ 

n 

-isinh^aXS^S^-S^-Si)} . (9) 

Using the Holstein-Primakoff creation and annihilation operators, expanding 
H(h) for large 5, keeping the 0(5 2 ) and 0(5) terms and performing a Bo- 

7 If T/A cut is small, the cut-off effects can be taken into account by including 4-derivative 
terms (with unknown couplings) in eq. (H). For p s 2> T ~ A cut , however, the systematic 
derivative expansion of chiral perturbation theory breaks down. 

8 The author is indebted to Uwe Wiese for explaining the proper way of introducing the 
chemical potential into the Heisenberg model. 
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goljubov transformation^ we obtained 

H(h)=Veo(h)+J2^k(h)blb k , [bk,bl] = 8 h ,k> , (10) 

where, in the infinite volume limitP| 

e (h) = -S{S+l)J{d+ l -{ha) 2 )+ l - ( -^u k (h), 



2 y ' J 2j (2tt) 2 
uj k (h) = 2JS[r k + {haf Sk \ 1 '\ (11) 
r k = — cos kj) /J(l + cosfcj) , 



- cos fci J^(l — cos fcj) + d , i,_j = l,2. 



We denoted by d and V the space dimension (d = 2 in our case) and volume, 
respectively. The momentum space integral runs over the Brillouin zone. The 
simple Hamilton operator in eq. (|l^) describes free magnon excitations and a 
zero-point energy. The magnon specie around k = (0, 0) picks up a mass ~ h 2 , 
the one around k = (-7T, 7r) remains massless since s k goes to zero there. It is an 
easy exercise to show that the contribution of the zero point energy to the free 
energy per unit slab area f(h) — /(0) is — h 2 pf w /(2T), where pf w is the SWE 
result for the spin-stiffness up to and including O(S). Actually, this should be 
so in any order of the SWE, due to the fact that p s is identical to the helicity 
modulus and h is an imaginary twist p^, p|. Adding the contribution from the 
magnon excitations we get: 



ffh\ ffn\ 1 h*P? W m 1 f d * k i l-exp(-^ fc (M) n9 , 

/ W - f = --h —— + — / ——^\n- —{ — . 12) 

T a 2 J (2tt) a 1 - exp(-^w fc (0)) 



We have to locate the 0(h 2 ) part of the integral in eq. ([12]). The result can be 
written in the following way: 

h 2 f o sw 1 hc sw 1 1 

m - /(») = - Y V + rJ»— - 5 - CWI • (13) 



where 



C{i) = J + In 8 + In 7 (14) 



2 

+ 27T7 2 



d 2 k s k rj k exp(-r; fc ) - 1 + exp(-7? fc ) 

(27r) 2 T) k T) k [l - eX P (-77fe)] |, )fc=7 (r fc )l/2 



9 These are standard manipulations in the literature on SWEp4]. 

10 Here and in the following we suppress powers higher than /i, since those terms in f(h) 
are are not universal pieces in the 2D cr-model. 
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using the notation 7 = 2JS/T. For T — * 0, exp(— (7(7)) behaves like 1 — 
1.05(1)7 -2 and eq. ( |l3| ) becomes identical with the spin-wave expanded form 
of eq. (7) in ref. [0| obtained from the effective theory. (There are additional 
cx (T/p s ) 2 corrections to the 0(T 2 ) term in eq. (||), of coarse.) This is an explicit 
confirmation of the effective action technique. 

The new term cx C(j) in eq. jl^ ) modifies the relation between the param- 
eters giving a = — C( , j)/(2tt) in eq. (0). Eq. (||) leads then to the result quoted 
in eq. (|). 

Let us close this paper with some remarks. It would be interesting to cal- 
culate the next order in the SWE which would test the full 2-loop result for 
consistency and add subleading cut-off corrections. This calculation seems to 
be feasible. 

In the large 5 limit the square lattice quantum Heisenberg model becomes 
identical to the standard lattice regularized 0(3) er-model for any correlation 
length. Due to the existence of powerful cluster MC techniques precise corre- 
lation length data are available from £ = 0(1) up to 0(200) (with finite size 
scaling techniques even beyond) |26| which can be compared with results in the 
quantum Heisenberg model from series expansion and MC for large 5. 

This paper is about cut-off effects in the quantum Heisenberg model. We 
should emphasize, this has nothing to do with the cut-off effects in the 2D 
non-linear er-model. We always considered large correlation lengths and so the 
cut-off effects in the er-model, which are suppressed as ~ £ -2 , are negligible. 
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Figure 1: The corrected theoretical prediction for the correlation length nor- 
malized by the asymptotic behavior (including the 0(T) corrections) in eq. (||) 
for spin 1/2, 1 and 5/2 (the upper, middle and lower solid lines, respectively). 
Correlation lengths larger than 10 lattice units are considered only. The MC 
data are from ref.Q [S = 1/2), ref.Jl5| (S = 1), the experimental points at 
S = 5/2 are taken from ref.||, while the circles come from series expansion fl7f. 
The dotted lines correspond to a semiclassical model p|| which performs well at 
small correlation lengths but does not have the correct asymptotic behavior. For 
Us, at S = 1/2 the values p s = 0.180 and c = 1.657|||, while at S = 1 and 
5/2 the SWE results were used. Similarly, 7 was connected to T/2irp s with 
the help of the SWE results. 
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